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Abstract: We present an accurate and efficient charge conserving electromagnetic 
particle-in-cell (EM-PIC) algorithm based on unstructured grids for the analysis of 
multipaction in coaxial cables. The present 2.5-D axisymmetric EM-PIC algorithm 
discretizes Maxwell-Vlasov equations and integrates the Furman-Pivi probabilistic 
model to take into account the secondary electron emission process. We adopt a local 
and explicit update for the Maxwell solver based on a finite-element time domain 
(FETD) scheme. The EM-PIC algorithm ensures charge conservation on 
unstructured grids from first principles by using Whitney forms. The present 
algorithm is used to investigate single- and double-sided multipactor effects in 
coaxial cables following multipaction susceptibility charts. In addition, we study 
surface treatments such as surface corrugations as possible multipaction mitigation 
strategies. This study is enabled by the use of unstructured grids that make it possible 
to capture fine geometric structures while preserving charge conservation. 
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U Introduction and motivation of our work 
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o Definition: 


e Electron resonance effect in presence of strong RF field, which causes discharge 
(emission) of electrons from conductor or dielectric surfaces. 

e Common in high power vacuum electronic devices or in some harsh environments 
(such as in spaceborne assets). 


o Practical scenarios of where multipactor may occur 
(undesired) 


e Harsh space weather conditions, vacuum electronic devices, high power microwave 
devices, RF satellite payloads, particle accelerators. 


o Practical scenarios of where multipactor can be 
exploited (desired) 


e Plasma displays, electron guns, dissipation of energy for receiver protection 





I. Introduction and Motivation of Our Work Multipactor 








o Snapshots of Multipaction Phenomena: (red dots represent superparticles) 


Conductor surface 


Conductor surface 





I. Introduction and Motivation of Our Work Multipactor 


Multipactor in Coaxial Cable 





e The radial electric field inside the coaxial cable 
is of the form, 


Outer Conductor 









Inner e 


e Asymmetric field generates Miller force which Iy PA a 


pushes the electrons towards outer conductor Single-sided | Double-sided 


o Single-Sided Multipaction 


e Due to the Miller force sometimes the secondary electrons emitted from the 
outer conductor fails to reach the inner conductor and travel back to the outer 
conductor, giving rise to the single-sided multipaction. 


ə Double-Sided Multipaction 


e Usual type of multipaction where the electrons emitted from one(inner/outer) 
conductor hits the other (outer/inner) conductor generating secondary electrons. 





R. Udiljak, D. Anderson, M. Lisak, V. Semenov, and J. Puech, “Multipactor in a coaxial transmission line. i. analytical study,” Physics of plasmas, vol. 14, no. 3, p. 033508, 2007. 
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Maxwell-Vlasov Equation 





o Multipactor effects are governed by, 
+ Maxwell-Vlasov system. 


+ Secondary Electron Emission process Lorant 


force 


o Maxwell-Vlasov equation 


describes the time evolution of collision-less plasmas as 






Spatial 
distribution 
of electrons 





Maxwell-Vlasov system 


dx dv 
Q a = Vp Q Mp = Fhorentz = dp [E (Xp) + Vp X B (Xp) 





E 
H = — +J 
V x a 


Particles are uncorrelated and collectively interact through long-range electromagnetic forces. 


Collisionless plasmas are hot (very fast) and dilute. 





II. Electromagnetic Particle-in-Cell (EM-PIC) Algorithm Maxwell-Vlasov Equation 


Superparticles and Meshing 








@ We use Particle-in-Cell (PIC) algorithm for solving Maxwell- 
Vlasov equation. 


Superparticle 


ə We are not interested in the microscopical 
behavior of each electrons, rather we are 
concerned about the overall behavior. So, we deal 
with superparticles which are coarse graining of 
continuous spatial distribution of electrons. S Electrons 


@ First, we discretize the space using unstructured triangular grids. 


e High geometric fidelity 

e Better grid dispersion 
performance. 
However, unstructured 
grid possess some 
unique problems for 


Structured Grid Unstructured Grid ) mam) EM-PIC (discussed 


later) 








II. Electromagnetic Particle-in-Cell (EM-PIC) Algorithm Superparticles and Meshing 


Basic Steps in EM-PIC 





Updates the particle 
(electron) position 
and velocity using 

Lorentz force 
equation and 
Newton's laws of 
motion 





o There are four fundamental steps in EM-PIC 



















Particle Pushe 


Xp Vp Yp 


Quantities not in mesh (continuous 






























Gather 


Interpolates fields from 
mesh to particle 
position 


Scatter 


/ Converts the motion and 
= position of particles to 
current and charges in 
mesh 










WHITN FORMS 








Updates the field 
values based on 
Maxwell’s 

equation 
















Quantities in mesh (discrete) 





II. Electromagnetic Particle-in-Cell (EM-PIC) Algorithm Basic Steps in EM-PIC 


Field Solver 











II. Electromagnetic Particle-in-Cell (EM-PIC) Algorithm Basic Steps in EM-PIC 


Differential Forms and Maxwell's Equation 





o Maxwell's equation in differential form is given by 
Electric field a Magnetic 
intensity "= _ flux density 
Ot 9 
ae Fur B) E Ot [te ) HI krear current 





Exterior derivative . 
density 
Primal | Dual 
Mesh Hodge operator Hodge operator Mesh 
taking into account taking into account 
the permeability of the permittivity of 
the medium the medium 
Differential Forms Vector Proxies 


€ = Enda + Eydy + Ezdz E = Ext + Eyj + E 
B = Badydz + Bydzdx + B¿dxdy «mu B = Byt + Byy + Bzz 
J = Sadydz + Jydzdx + J¿dxdy I = Jal + Jyy + J¿2 








II. Electromagnetic Particle-in-Cell (EM-PIC) Algorithm Field Solver 


Field Expansion and Interpolants 





in primal mesh 


Nı 
E (r,t) ~ XC E; (t) W9 (r) Primal Mesh 


I= 


B (r,t) x SB, (t) WO (r) 
k=l 





N2 
J, (r,t) ~ dates (t) wi? (r) Bt? = [B]? — At [Deun] : [E]" 


D (t) WA (r) * el ° E+ = [x el * Ej” + At[D Deurt| | By" +2 4 nts 


Leapfrog time integration 





Discrete Hodge matrices. 
(constitutive relations and all 
metric info.) 








II. Electromagnetic Particle-in-Cell (EM-PIC) Algorithm Field Solver 


Basic Steps in EM-PIC 








II. Electromagnetic Particle-in-Cell (EM-PIC) Algorithm 


Field Interpolation via Whitney Forms 





ə Interpolation of electric field at particle's position from the mesh edges 
via Whitney 1-forms: 





In a similar fashion, Bẹ can be 
evaluated with Whitney 2-forms. 





Total field 





II. Electromagnetic Particle-in-Cell (EM-PIC) Algorithm Gather 


Basic Steps in EM-PIC 


Particle 
Pusher 








II. Electromagnetic Particle-in-Cell (EM-PIC) Algorithm Particle Pusher 


Particle-Pusher (PP) 





@ Several relativistic PPs : (1) relativistic Boris, (2) Vay, (3) Higuera-Cary 





II. Electromagnetic Particle-in-Cell (EM-PIC) Algorithm Particle Pusher 


Basic Steps in EM-PIC 


II. Electromagnetic Particle-in-Cell (EM-PIC) Algorithm Scatter 








Evaluation of Grid Charge and Current 





o Charge is assigned to nodes while current is assigned to edges 
again based on Whitney O- and 1-forms. 








o Issue of Charge 
Conservation 

o Needs correction 
terms to enforce 

charge conservation 






Structured Grid = Unstructured Grid 


@ Our proposed algorithm is able to maintain exact charge conservation 
as it uses Whitney forms as interpolatory functions. 


o Exact charge conservation is achieved from first principle without any 
need of correction terms. 





H. Moon, F. L. Teixeira, and Y. A. Omelchenko, “Exact charge-conserving scatter—gather algorithm for particle-in-cell simulations on unstructured grids: A geometric 
perspective,” Computer Physics Communications, vol. 194, pp. 43-53, 2015. 


II. Electromagnetic Particle-in-Cell (EM-PIC) Algorithm 





Scatter 


Evaluation of Grid Charge and Current 





ə Charge is assigned to nodes while current is assigned to edges 
again based on Whitney 0- and 1-forms. 


Grid charge: Q) = qC> Ww (x p) = = GpCp Aj 


mm LAA — pe 


n+ O A 
Grid current: E = we f wo (x) - dl 


during At 


Barycentric coordinate system 
An, 
An, = vr a A 


Any = An, {Ar 


3 
T 2 Am =4 


== (Ars ) Nae ee) 





Point weighting Line weighting 
on node on edge 





exact charge conservation 





II. Electromagnetic Particle-in-Cell (EM-PIC) Algorithm Scatter 
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Body of Revolution 










Axisymmetric 


cine, > m=0,+1, +2.... 


Fourier mode 
oe A 0 expansion 


Cross section of coaxial cable 












No azimuthal 
variation 





Due to cylindrical symmetry of 
multipaction phenomena, our 
problem reduces from 3-D to 2-D. 


pz plane is discretized in unstructured 


triangular mesh 


HI. Adaptation to Axisymmetric Case Body of Revolution 


Artificial Constitutive Parameters 





o We can correlate the TM, case (E,,E,,B,¿) with TE, case (E, ,E,,B,). 
o TE,case (E, ,E,,B.): 


= *,-18 =w*Brdz 


D = xE = e(Endydz + Eydxdz) 


p> l Cross section of coaxial cable 
o TM, case (E,,£,,B4): H 





H = *,-18 =(u PB gdp oe Artificial constitutive 
D= tet Eos + Esdpdo) parameters 


/ — ep 


€ 





H =x mB = W7 Bado 
D = xE = e'(Epdodz + Ezdpdo) 
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Boundary Conditions for Charged Particles 





Macroscopic components in the model 
Backscattered (BS) : elastic 


9 Secondary electron emission (SEE) | Re-diffused (RD) : partially elastic 


True-secondary (TS) : inelastic 





= Furman - Pivi probabilistic model based 
on experimental data for various materials. 


= (i) Secondary electron yield (SEY), 
(11) emitted-energy spectrum 


metallic surface 


= Monte-Carlo method depending on 
incident energy and angle 


e M. Furman and M. Pivi, “Probabilistic model for the simulation of secondary electron emission,” Physical Review Special Topics-Accelerators and Beams, vol. 5, no. 12, p. 


124404, 2002. 


e D.-Y. Na and F. L. Teixeira, “Analysis of multipactor effects by a particle-in-cell algorithm integrated with secondary electron emission model on irregular grids,” IEEE 


Transactions on Plasma Science,vol. 47, no. 2, pp. 1269-1278, 2019. 
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Problem Parameters 
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line current 
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| 1 Longitudinal 


cross-section of 


coaxial cable 1 GHz 


Outer conductor 







Inner conductor 


Inner conductor 


800 V 
Outer conductor 


Single-sided 
multipactor 


Unstructured 
triangular mesh 





Double-sided multipactor 


Fxcitation voltages are chosen based on multipaction susceptibility chart. 





V. EM-PIC Simulations of Multipaction in Coaxial Cable Problem Parameters 
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Simulation Results 
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Red dots represent superparticles 





V. EM-PIC Simulations of Multipaction in Coaxial Cable Simulation Results 


Simulation Results 
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IV. PIC SIMULATIONS OF MULTIPACTOR Simulation Results 


Conclusion 


o The integration of a charge- and energy-conserving finite- 
element-based PIC algorithm was implemented on unstructured 
meshes with Furman-Pivi probabilistic model describing 
secondary electron emission processes on metallic surfaces. 


o The algorithm was employed to analyze multipaction (both single 
and double-sided) in coaxial cables. 


@ The multipactor saturation process for single-sided case was 
examined by quantifying the output power loss and harmonic 
generation arising from acceleration phase mismatching and self- 
field counterbalance effects. 


@ The Proposed SEE-PIC algorithm enabled mesh refinement and was 
better suited to model to complex geometries. The algorithm was 
employed to observe multipaction mitigation in grooved coaxial 
cable. 
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